Stop searches in 2012 



IPPP/12/30 
DCPT/12/60 



Tilman Plehn, 1 Michael Spannowsky, 2 and Michihisa Takeuchi 1 

1 Institut fur Theoretische Physik, Universitat Heidelberg, Germany 
2 IPPP, Department of Physics, Durham University, United Kingdom 



For this year's 8 TeV run of the LHC we lay out different strategies to search for scalar top pairs. 
We show results for the hadronic and for the semi-leptonic channels based on hadronic top tagging. 
For the di-lepton channel we illustrate the impact of transverse mass variables. Each of our signal- 
to-background ratios ranges around unity for a stop mass around 400 GeV. The combined signal 
significances show that dedicated stop searches are becoming sensitive over a non-negligible part of 
parameter space. 
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I. INTRODUCTION 



The primary goal of the LHC is to unveil the mechanism which breaks electroweak symmetry. The minimal 
solution, realized in the Standard Model (SM), predicts the existence of a fundamental scalar field, the Higgs 
boson. The perturbative instability of scalar masses is one of the motivations for the existence of new physics at 
the electroweak scale [TJ . One of the proposed solutions to this problem is supersymmetry at the TeV scale [2] • 

The leading contribution to the quadratic divergence of the scalar Higgs mass in the Standard Model arises 
from the top quark with its large Yukawa coupling. Its supersymmetric partner is the scalar top quark. 
It can be pair produced at hadron colliders just through its QCD couplings, i.e. without any mediation by 
additional supersymmetric particles. Its production cross section at the LHC therefore only depends on the stop 
mass [3] . The decay channels of the stop depend almost entirely on the supersymmetric mass spectrum [H [5] : 
if kinematically allowed, it will decay through its strong coupling into a top quark and a gluino. Below this 
threshold there exist two weak decay channels, a charged current decay to bottom-chargino and a neutral 
current decay 

h -> txi , (i) 

where we assume the lightest neutralino to be the dark matter agent and hence appear as missing transverse 
momentum at hadron colliders. This neutral current decay we target in this analysis. For even smaller stop 
masses a loop-induced decay to a charm quark and the lightest neutralino dominates. 

At the Tevatron searches for such loop-induced decays as well as the bottom-chargino signature have lead to 
moderate limits on the stop mass [51 [7]- In 2011 the LHC delivered approximately 5 fb _1 of integrated luminos- 
ity to ATLAS and CMS. This data set allowed both collaborations to significantly constrain the squark-gluino 
mass plane through the search for jets plus missing energy [9]. 1 Direct stop searches at the LHC are 
notorious, mostly because of the small signal cross section and the overwhelming background from top pair 
production [llj . Therefore, all limits on top partner searches from 2011 analyses are derived from supersym- 
metric model assumptions [12j : either, the gluino has a sizable branching ratio into 3rd generation squarks |13j 
or sbottom limits get re-interpreted in terms of the stop mass based on the SU(2) symmetry of their left-handed 
modes. More specifically, assuming gauge mediated supersymmetry breaking with the NLSP decay xl ~ * ZG 
we can exclude masses below < 330 GeV for m^o = 190 |14) . 

For 14 TeV collider energy we have shown that with the help of top taggers jT5J HFJ we can reconstruct 
top quarks from stop decays and extract stop pair production from top backgrounds. This holds for purely 
hadronic decays of the two top quarks from the stop pair [TTJ [17] as well for the semi-leptonic channel [18] . 
The main benefits of the new top reconstruction methods in these analyses is that they automatically resolve 
any combinatorics of the top decay jets and fully reconstruct the top 4-momenta and angular correlations [T§] 
The associated analysis is not any more complicated than a search for slepton or sbottom pairs, including the 
application of tot2 for a stop mass measurement [20] . 

In this paper we test a variety of search strategies for light stops (350 < < 700 GeV) at 8 TeV. We see 
that a straightforward adaption of the 14 TeV search strategies is challenging, due to the small signal rate. 
Instead, we can optimize the signal efficiency by looser requirements on the reconstructed stop decays while 
benefiting from the also significantly reduced top pair background. This change with respect to the 14 TeV 
analyses [TTJ [T71 [TH] reduces the potential of reconstructing masses and model parameters, but it will allow us 
to probe a significant range of stop masses in 2012 [2T|. 

We organize this brief overview in the following way: In Sec. [TTJ we briefly review the properties of the stop 
signal and of the SM backgrounds at 8 TeV. This includes production cross sections and px distributions. In 
Sec. |III| we discuss our results for hadronic, semi-leptonic and leptonic top decays. In the latter case we find 
that a simple transverse mass variable removes almost the entire background. For each of these channels we 



quote signal-to-background ratios and signal significances. A summary follows in Sec. IV 



1 Note that replacing the squark-gluino mass plane for example by an nio-inyj plane in the CMSSM implies significant unwanted 
and unnecessary model assumptions [101 . 
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II. STOP PRODUCTION AT 8 TEV 



Although the LHC production cross section for stop pairs is much smaller at 8 TeV than at 14 TeV there 
is a fair chance to disentangle this signal from the also significantly smaller SM backgrounds. For the rest of 
the paper we omit the index '1' for the lighter of the two stops. Stop pair production always refers to the pair 
production of the lighter of the two stop mass eigenstates. In Fig. [I] we show the NLO signal cross sections for 
y^s = 8 TeV and 14 TeV as a function of the stop mass [22]. It drops by roughly a factor 1/6 from 14 TeV to 
8 TeV collider energy. This reduction is partly compensated for by the tt+jets cross section, which is reduced 
by a factor 1 /4. We consistently normalize our tt cross section to the approximated NNLO result of 234 pb at 
8 TeV and 918 pb at 14 TeV [23]. 

In our 14 TeV analyses we have shown that relying on boosted hadronic tops we can achieve S/B ~ 1 and 
S/\/~B > 5 in the fully hadronic mode [11] and S/B ~ 2 and S/^fB > 5 in the semi-leptonic mode [18] for 
10 fb -1 . For the 2012 run at 8 TeV the envisioned integrated luminosity also ranges around 10 — 20 fb^ 1 . With 
the small production rates at 8 TeV we need to adapt our search strategies to retain enough signal events to 
achieve a good signal significance S/\[B for 10 fb -1 . 

To decide on how to improve the analyses it is instructive to study kinematic correlations between the two 
tops in the signal process 

pp -> tt* -> (tx?) (ffi) (2) 

and the tt background. Fig. [2] shows the top versus the anti-top transverse momenta for the signal and 
background, assuming mi — 400 GeV. For the it background we find a strong correlation between the top 
and the anti-top, while for tops from stop decays there is almost no correlation. The initially back-to-back 
configuration of the stops is significantly distorted by the \i LSP momenta. 

A sizable (transverse) top momentum is the key to top tagging. In Tab. [I] we show the predicted rates for 
at least one top with px.t > 200 GeV. This is the value we need to reliably apply subjet techniques and fully 
reconstruct hadronic tops |17) . For the top pair background we include up to two hard jets when we simulate 
the distributions. We see that these additional jets lead to a significant fraction of events where only one of 
the two tops drops below 200 GeV. The signal with its comparably uncorrelated tops clearly prefers only one 
strongly boosted top quark. 

The last column in Tab. [I] shows the ratios of rates for 400 GeV stops divided by the tt background rate. 
Indeed, we gain a factor 2.7 in S/B by focusing on events with only one top with pr jt > 200 GeV. This means 
that focusing on events with asymmetric p?t by only requiring one top tag should be the strategy for the 2012 
run at 8 TeV. In the following subsections we study four statistically independent samples or analyses: 

• two hadronic boosted tops 

• one hadronic boosted top and one hadronic un-boosted top 

• one hadronic boosted top and one leptonic top 

• two leptonic tops 
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Figure 1: Production cross section for (s)top pairs as a function of the mass m t - assuming \/S ■ 



14 and 8 TeV. 
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Figure 2: Top vs. anti-top transverse momenta for 400 GeV stop pairs (left) and the tt background (center) at 8 TeV. 
The right panel shows the ratio tt* /tt. 



To generate the signal sample we rely on Herwig++ [33] and assume stop masses of = 350, 400, 450, 
500, 600, and 700 GeV. We normalize their production cross section to the Prospino results at next-to-leading 
order [3J. They range from 0.8 pb (mj = 350 GeV) to 0.1 pb (500 GeV) as shown in Fig.JT] Resummation slightly 
increases this rate further [35]. For all stop masses we choose the same lightest neutralino mass m^o = 100 GeV. 

Our Standard Model backgrounds are tt+jets, QCD jets, VF+jets, Z+jets, and ttZ. We use Alp- 
GEN+Pythia [35] [33] to generate the corresponding samples. For all processes except for ttZ production 
wc use MLM matching [37] to simulate additional hard radiation. We match up to ti+2 jets, W+3 jets, 
Z+4 jets, and 3 — 5 jets for the QCD sample. 

The leading tt+jets background sample we normalize to the approximate NNLO rate of 234 pb [23]. For the 
subleading background channels we use the leading order normalization. The ttZ cross section at LO yields 
21.5 fb, based on Alpgen including the Z —> vv branching ratio. Since the tiZ rate only becomes comparable 
to the stop rate for > 600 GeV we neglect this irreducible ttZ background. We have checked that for all 
processes considered this does not affect the quoted results. 

Our analysis is based on a simple calorimeter simulation with granularity of 0.1 x 0.1 in (t],4>). We sum 
the four momentum of all particles in each cell and rescale the resulting three-momentum such as to make 
the cells massless. The calorimeter cells are later on used as (fat)-jet constituents. Throughout this work 
we use the Cambridge/ Aachen (C/A) algorithm [3S] with R — 1.5, as implemented in FastJet [35]. The 
resulting fat jets are then used as input for the HEPTopTagger. Preliminary ATLAS analysis show that the 
HEPTopTagger results are only very mildly affected by detector effects, underlying event, or pile-up [3D]. 
For regular QCD jets use the same C/A algorithm with R = 0.5. When analyzing leptonic or scmilcptonic top 
decays we require the leptons to be hard and isolated: pt,i > 20 GeV and i^T,had < 0.1 Ere within R < 0.2 
around the lepton. 
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Table I: Signal and background cross sections [fb] for different stop masses. We assume BR(i — > txi) = 1. 
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III. STOP RECONSTRUCTION AT 8 TEV 

As mentioned in the last section, we present four distinct signatures of stop pairs decaying to a pair of top 
quarks plus missing energy. The different analyses are chosen such that they are statistically independent and 
can be combined, if required. 

A. Two top tags 

As a first attempt to apply the successful 14 TeV strategy to the 2012 run, we analyze fully hadronic stop 
pairs, 

pp -+ it* -> (ix°) (ixl) -> (bjjjfi) (bjjxl) ■ (3) 

One of the advantages of the fully hadronic mode is that, in principle, the top tagger fully reconstructs both 
top momenta. This means that the fully hadronic channel is particularly well suited for detailed studies of a 
top partner signal. For 14 TeV we indeed expect signal-to-background ratios around S/B ~ 1 [11]. In this 
section we follow the same analysis steps at 8 TeV. 

For the fully hadronic double tag mode we veto isolated leptons and require at least two fat jets with 

Pt.j > 150/150 GeV and i> T > 100 GeV , (4) 

where both numbers are slightly reduced compared to the 14 TeV case [TT] . 

In Tab. [n] we show the numbers of events after each cut. For illustration purposes we compute S/y/B for 
10 fb" 1 after each step, taking into account only statistical errors. 

We first require missing transverse momentum, just like for almost all supersymmetry analyses. To simulate 
the amount of fake missing energy in the pure QCD jets background we would need a detailed detector simula- 
tion. Instead, we adopt a conservative efficiency of 1% for QCD events to pass this cut. With this assumption 
QCD jets are our dominant background, which essentially forces us to apply a 6-tag inside the first tagged top 
later on. 

This analysis is based on two top tags, using the HEPTopTagger algorithm. Already at this second stage 
the IV+jets and Z+jets backgrounds become negligible and we are left with it and QCD jets. On the other 
hand, all rates are already down to a femto-barn level. As mentioned above, we need an additional tagged 
6-subjet inside the first tagged top, to reduces the QCD background to a negligible level. It has been shown 
that ^-tagging is not a problem inside a fat jet, so we apply a 6-tagging efficiency of 50% for 6-quark and a 1% 
fake rate for light quarks and gluons [31] . In the final step the ti background can be reduced by requiring 

m T2 > 250 GeV . (5) 

This way we achieve S/B = 0.79 and S/yB = 1.5 for an integrated luminosity of 10 fb" 1 . The fact that the 
significance is reduced by a factor 1/3 as compared to the 14 TeV run is explained by the factor 1/6 in signal 
cross sections and the softer px,t spectrum. 
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Table II: Analysis flow for the two-top analysis. All numbers are given in fb. The symbol "— " denotes less than 0.01 fb. 
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To study ways out of this rate limitation we can turn to the top tagging efficiencies. After requiring $t > 
100 GeV we show these efficiencies in Tab. [TTJ For the first (mis-)tag it ranges around 30% for the signal, 20% 
for ti, and 1% to 2% for QCD jets, IV+jets, and Z+]ets. For the second tag the efficiencies are around 10% 
for the signal, 2% for ti, while QCD jets, W+jets, and Z+jets remain at 1 to 2%. The reason for this small ti 
efficiency is that in events with large missing momentum one of the tops has to decay leptonically. The second 
tag then is a fake-top from the remaining b jet combined with hard QCD radiation. Hence the second top tag 
is very helpful against the leading top pair background. Unfortunately due to the already small signal rate, the 
second top tag does not increase the sensitivity significantly. 



B. One top tag and one bottom tag 

To improve the fully hadronic analysis presented in the last section we propose a search for one boosted top 
and one &-tag in the recoiling softer top decay jets. As a starting point we apply a lepton veto and require 
exactly one fat jet together with missing transverse momentum, 

pr,j > 150 GeV and fa > 100 GeV . (6) 

One subjet inside the tagged top has to be 6-tagged. In addition, we require a continuum 6-tag which cannot 



be a constituent of the tagged top. In Tab. [HI] we see that after these two 6-tags all backgrounds except for 
top pair production are negligible. 

To reduce the still overwhelming ti background we construct a specific transverse mass variable from the 
general form 

'vis 



mr(priB,pr) = \ mj is + 2\fa\(E T y ls -p TjVis cos0) , (7) 



where <f> is the angle between the transverse visible and missing momenta. In events with an isolated lepton 
and missing momentum this variable is commonly used to reject leptonic W decays because mj- is bounded 
from above. 

The main background in the no-lepton mode with large missing momentum comes from ti events where one 
top decays through a tau lepton. In these events large missing momentum can be induced by the neutrinos 
from the W and r decays. To separate these events from the signal with its two neutralinos we construct the 
transverse mass with a 6-jet instead of the lepton, and require 

m$ = m T (pb,pr) > 200 GeV. (8) 

The left panel of Fig. [3] shows the normalized m^' distributions for the signal and the ti background. As 
expected, there is an endpoint at m* for the background. As a cross check the thin dashed line shows the 
mSP distribution from the missing momentum and the bottom momentum from a leptonic top decay in semi- 
leptonic ti events. This parton level distribution shows good agreement with the background distribution from 
the hadronic final state. We also tested a similar m^\, but the remaining number of signal events turns out to 
be too small. 
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Table III: Analysis flow for one top tag and one 6-tag. All numbers are given in fb. The symbol "— " denotes less than 
0.01 fb. In the last line we illustrate the potential of a 100% efficient tau veto. 
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Figure 3: Left: 



distribution for the hadronic signal (solid) and the tt background (dashed). The thin dashed line 



shows tt semi-leptonic events at parton level. Right: mr distribution for the signal, for the semi-leptonic top pairs and 
for di-leptonic top pairs. The distributions are based on events which have exactly one lepton (e, fx) in the final state. 



After this final 



cut all signal and background rates shown in Tab. 



Ill are again at the fb level. For 

a 400 GeV stop we find S/B = 0.73 and S/y/B = 2.6 with 10 fb -1 of data. Because we know the origin of 
the remaining tt events we could further improve the results by rejecting tau leptons. Just to illustrate the 



potential of such a requirement, we assume a 100% efficiency for tau rejection in the last line of Tab. Ill 



C. One top tag and one lepton 

Following the previous section, events with one boosted and one non-boosted top are well suited to extract 
stop pairs from Standard Model backgrounds. An obvious question then becomes what happens if the softer 
of the two tops decays leptonically, 

pp it* (tx?) (txS) {U + vxD (bjjxl) + (bjjxl) {bTuxl). (9) 
This time we require one isolated lepton, a sizable amount of missing energy and one fat jet with a top tag. In 



Table IV we see that the tt background is still overwhelming. The reason is a significant fraction of semi-leptonic 
tt events passing these cuts. 

The transverse mass, Eq.([7]), has an upper kinematic endpoint for events where the missing energy comes 
from leptonic W decays. To efficiently reject leptonic top pair events as well as any kind of W events we 
require [32] 

m T > 150 GeV . (10) 
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Table IV: Analysis flow for one top tag and one lepton. All numbers given in fb. The symbol "— " denotes less than 
0.01 fb. 



8 




Figure 4: Normalized m" 2 distributions for the stop signal and the tt background. The light dotted line for the tt 
background includes a Gaussian smearing of the missing energy. 



Surprisingly, a non- negligible number of tt events survives this cut. Similarly to the last section, purely leptonic 
top pairs can fake a top tag from the b jets and additional QCD radiation. The right panel of Fig. [3] shows the 
normalized tut distributions for the signal, for semi-leptonic tt events, and for purely leptonic ti events. 

After imposing all the above cuts we arrive at a promising signal-to-background ratio of S/B = 1 and a 
significance of S/yB = 6.5 for 10 fb _1 of 8 TeV running. Unlike the hadronic channels this analysis does not 
require a b inside or outside the top tag. However, if we are willing to pay the price in available rate we can 
apply the usual 6-tag among the top tag constituents. 



D. Two leptons 



Until now, all our stop pair analyses involve one boosted hadronic top decay which we identify and reconstruct 
using a top tagger. If we loosen our requirements on event reconstruction we can of course search for top pairs 
in purely leptonic top pairs, 

pp -> tt* -> (tx?) (*X?) "> (bt + vXi) (bruxl) . (11) 

This di-lepton channel turns out to have the largest reach in finding or ruling out anomalies in the top sector. 
On the other hand, in the absence of any reconstructed mass a deviation from the Standard Model cannot 
confirm the existence of top partners. Therefore, we consider the di-lepton channel a very powerful tool to 
confirm and statistically support any anomaly found in one of the hadronic channels. 

The previous sections show that the transverse mass vtlt with a lepton or a fo-jet momentum efficiently reduces 
semi-leptonic tt backgrounds. For events with two sources of missing energy a better-suited variable is |20j 

m" 2 = min [max {m T (pe x , $i ), m T (pt 2 , fa))} ■ (12) 
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m? 2 > 100 GeV 
m^ 2 > 150 GeV 


6.05 4.30 2.70 1.65 0.56 0.20 
0.81 1.21 1.06 0.81 0.34 0.14 


0.65 (0.79) 0.09 
0.00 (0.03) 0.02 


5.8 (4.9) 15.8 (14.5) 
n.a. n.a. 



Table V: Analysis flow for the di-lepton mode. All numbers are given in fb. The ttZ numbers are shown including the 
decay Z — > vv. The number in parentheses include a smeared transverse momentum measurement. 
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In this definition the two transverse mass values, Eq.([7]), are computed with one of the two observed leptons. The 
missing momenta in the two transverse mass values are a hypothetical split of the observed missing transverse 
momentum into two parts. The mass of the unobserved particles we assume to vanish. 

To enhance the di-lepton mode we need to identify additional sources of missing energy, as compared to 
tt or W + W~ events with their W decay neutrinos. The particles contributing to m,T2 from tt events come 
from the W + W~ subsystem, so for top quarks close to threshold the upper endpoint is tuti < mnw O 151] . 
In contrast, for the stop signal such an endpoint does not exist. The corresponding 771^2 distributions we 
show in Fig. [4 | Th e leading detector effect is the smeared missing momentum according to a Gaussian with 
a(fi) = a ■ ^Y^E T ,a = 0.53 - 0.57 @S]. The thin dotted line in Fig. [I] includes this smearing and shows a 
slightly enhanced tail. 

Our purely leptonic event selection starts with two isolated leptons and fir > 100 GeV. Tab. [V] shows that 
the tt background becomes essentially negligible after we require 

m^ 2 > 100 GeV. (13) 

Similarly, the total ttZ rate is significantly smaller than the signal. Including a smeared missing energy 
measurement increased the tt background from 0.65 to 0.79 fb. For = 400 GeV, we expect more than 40 
signal events, giving S/B = 5.8 and S/v~B = 15.8 with 10 fb -1 . A slightly harder cut m" 2 > 150 GeV removes 
essentially all the SM backgrounds. What is most impressive is that this analysis completely ignores any jet 
which might come with the £ + £~fir system! 

IV. SUMMARY 

Reasonably light top partners are necessary to solve the hierarchy problem. Therefore, searches for stops or 
other top partners are of paramount interest to LHC physics. In 2012 the LHC will gather at least O(10) fb -1 
of data at 8 TeV. For four independent search channels we show how 2012 data will start to either find or 
exclude light top partners, decaying to top quarks and missing energy. 

In the fully hadronic mode we study two strategies: tagging either one or two hadronic tops we find S/B ~ 1 
for a stop mass of 400 GeV. Unfortunately, the statistical significance is rather modest, S/V~B = 1.5 (two tags) 
and S/VB = 3.0 (one tag). 

Searches for semi-leptonic or fully-leptonic top pairs are more promising. In the semi-leptonic mode we tag 
one top recoiling against an isolated lepton. After cutting on mr we find S/B = 2.1 and S/VB = 5.4. In 
the di-lepton mode a cut on m" 2 rejects almost all Standard Model backgrounds. This gives us a striking- 
sensitivity of S/B = 5.8 and S/^fB = 15.8. 

Obviously, the fully leptonic mode is unlikely to conclusively reconstruct and confirm a top partner. However, 
the combination with the statistically less significant hadronic modes should allow us to establish a top partner 
signal in 2012. 
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